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ABSTRACT 

Galaxy clusters contain a large population of low mass dwarf elliptical galaxies whose exact 
origin is unclear: their colors, structural properties and kinematics differ substantially from 
those of dwarf irregulars in the field. We use the Illustris cosmological simulation to study 
differences in the assembly histories of dwarf galaxies (3 x 10® < M*/Mo < 10i°) accord¬ 
ing to their environment. We find that cluster dwarfs achieve their maximum total and stellar 
mass on average ^ 8 and ^ 4.5 Gyr ago (or redshifts z — 1.0 and z = 0.4, respectively), 
around the time of infall into the clusters. In contrast, field dwarfs not subjected to environ¬ 
mental stripping reach their maximum mass at z = 0. These different assembly trajectories 
naturally produce a color bimodality, with blue isolated dwarfs and redder cluster dwarfs ex¬ 
hibiting negligible star-formation today. The cessation of star-formation happens over median 
times 3.5-5 Gyr depending on stellar mass, and shows a large scatter (~ 1-8 Gyr), with the 
lower values associated with starburst events that occur at infall through the virial radius or 
pericentric passages. We argue that such starbursts together with the early assembly of cluster 
dwarfs can provide a natural explanation for the higher specific frequency of globular clusters 
(GCs) in cluster dwarfs, as found observationally. We present a simple model for the forma¬ 
tion and stripping of GCs that supports this interpretation. The origin of dwarf ellipticals in 
clusters is, therefore, consistent with an environmentally-driven evolution of field dwarf irreg¬ 
ulars. However, the z = 0 field analogs of cluster dwarf progenitors have today stellar masses 
a factor ~ 3 larger -a difference arising from the early truncation of star formation in cluster 
dwarfs. 

Key words: galaxies: dwarfs - galaxies: evolution - galaxies: interactions - galaxies: star 
clusters: general - methods: simulations 


1 INTRODUCTION 

By number, dwarf elliptical galaxies (dEs) dominate the population 
of dwarfs in cluster environments. They are, however, almost ab¬ 
sent in the field, where the prevalent morphological type is gas-rich 
dwarf irregulars (dirrs) (Binggeli et al. 1990; Grebel 1999; Geha 
et al. 2012). This difference immediately raises questions about the 
relationship between dEs and dirrs. Are these two types of galax¬ 
ies intrinsically different objects? Or, is there an evolutionary link 
between these two dwarf populations? 

There is no strict definition for dE galaxies, but they are rather 
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characterized by a combination of attributes such as having low 
mass, being of low surface brightness, and exhibiting red colors and 
low levels of star formation. Moreover, dEs have quite a heteroge¬ 
neous morphological mix (Lisker 2009; Janz et al. 2012) and also 
show wide kinematic variations (Geha et al. 2002; van Zee et al. 
2004a; Lisker2006; Rys et al. 2013; Toloba et al. 2014a,b). These 
observations hint at a complex formation scenario where more than 
one factor determines their final fate. Although the definition is not 
specific, there is clear consensus that dEs are less-disky and less 
rotationally-supported than dwarf irregulars in the field. 

Differences between the two populations could arise from ei¬ 
ther internal or external mechanisms. It appears that dirrs and dEs 
follow similar scaling relations, so they could be the same type of 


© 0000 RAS 


2 Mistani et al. 


object but in which internal processes, like star-formation and feed¬ 
back, caused some to retain more gas than others (Dekel & Silk 
1986; de Rijcke et al. 2005). While appealing, such a scenario can¬ 
not explain the strong correlation observed with the environment. 
Alternatively, dEs may have been dirrs in the past that were trans¬ 
formed into dEs later, due to environmental effects within their host 
clusters. Such a transformation could involve some degree of ram- 
pressure stripping and gas starvation to redden their colors (Con- 
selice et al. 2003, van Zee et al. 2004b, Boselli et al. 2008), and 
tidal effects to heat them kinematically, altering their morpholo¬ 
gies (Mayer et al. 2001; Gnedin 2003; Lisker 2009; D’Onghia et al. 
2009; Smith et al. 2010; Lelli et al. 2014). A different and intriguing 
perspective for the influence of environment on dwarfs was argued 
by Sabatini et al. (2005), who noted that interactions with the clus¬ 
ter medium could temporarily enhance star formation in infalling 
dirrs, leading to rapid consumption of the gas and a transformation 
into dEs. Finally, low mass dEs could also be the stripped remnants 
of more massive galaxies that fell into the cluster long ago (Con- 
selice et al. 2003). 

Environmentally driven transformation models for dEs, al¬ 
though attractive, are not free from challenges. One of the main 
problems with this scenario is the observed difference in the fre¬ 
quency of globular clusters (GCs) between dwarfs in clusters and 
in the field. If dEs and dirrs are evolutionarily linked, dEs should 
contain similar or even fewer GCs than disky galaxies in the field 
because of tidal stripping within the cluster host. Instead, dEs in 
clusters tend to show an increased GC specific frequency over dirrs 
at similar masses (Miller et al. 1998, Miller & Lotz 2007, Jordan 
et al 2007b, Peng et al. 2008, Sanchez-Janssen & Aguerri 2012). 
This discrepancy appears accentuated for nucleated dEs (Miller & 
Lotz 2007) and dwarfs in the inner regions of clusters (Jordan et al. 
2007b, Peng et al. 2008). 

With their capability of being able to track the evolution of 
galaxies, numerical simulations could, in principle, resolve these 
issues. In practice, this is difficult because GC formation is still 
poorly understood cosmologically (however see an original at¬ 
tempt by Prieto & Gnedin 2008), and the dynamic range needed 
to self-consistently evolve dwarf galaxies with their GC systems 
in a galaxy cluster environment is currently prohibitive. Here, we 
propose a hybrid approach. Hydrodynamical simulations of a large 
cosmological volume are performed to enable us to sample a set 
of galaxy clusters and, simultaneously, to follow the assembly and 
evolution of dwarfs in clusters and in the field. In post-processing, 
we then apply a simplistic model of GC formation together with 
a particle-tagging technique to roughly describe the formation and 
stripping of GCs in cluster dwarfs. Our effort builds on the earlier 
work by Peng et al. (2008) based on semi-analytical catalogues and 
incorporates the effect of tidal stripping of GCs not addressed in 
their work. 

This paper is organized as follows. We introduce the simu¬ 
lations and selection criteria for our samples in Sec. 2. The as¬ 
sembly of dwarfs according to their environment is investigated in 
Sec. 3 and the details of their star formation histories is presented 
in Sec. 4. We develop and apply our method for tracking GCs in 
Sec. 5, addressing the expected GC frequencies for cluster and field 
dwarfs. We summarize our main results in Sec. 6. 


2 NUMERICAL SIMULATIONS 


We employ a set of model galaxies taken from the Illustris 
Simulation (Vogelsberger et al. 2014a,b; Genel et al. 2014). Illustris 
consists of a cosmological box 106.5 Mpc on a side evolved with 
hydrodynamics using the code AREPO (Springel 2010). The simu¬ 
lation adopts a WMAP9-consistent cosmology and includes treat¬ 
ments of gravity together with the most relevant physical processes 
driving galaxy evolution, including gas cooling/heating (Katz et al. 
1996; Faucher-Giguere et al. 2009), star-formation (Springel & 
Hernquist 2003), metal enrichment, and AGN feedback (Springel 
et al. 2005; Di Matteo et al. 2005; Sijacki et al. 2007). While a 
detailed description of the model can be found elsewhere (Vogels¬ 
berger et al. 2013), for completeness we briefly summarize the rel¬ 
evant aspects of the physics included in our simulations. 

Gravity is computed using a Tree-PM approach in which long- 
range forces are calculated on a grid while short-range interactions 
are obtained by expanding the potential into multipoles (Springel 
2005). The Euler equations governing the hydrodynamics of the 
gas are solved on an unstructured mesh that moves with the flow, 
yielding a continuously adaptive method to achieve high resolu¬ 
tion where needed and simultaneously moderating the impact of 
numerical diffusion (Vogelsberger et al. 2012; Sijacki et al. 2012; 
Genel et al. 2013; Nelson et al. 2013). In this scheme, the gas is 
distributed among “cells ”, which can be heated or cooled accord¬ 
ing to the local density and metallicity. Dense gas becomes eligible 
for star-formation above a density threshold nn = 0.13 cm~®, at 
which point star particles are spawned probabilistically according 
to the local dynamical time. 

We follow the evolution of stars, including mass and metal 
return from stellar winds and the explosion of massive stars (su¬ 
pernova type I and II; hereafter, SNI and SNII, respectively). In 
the case of SNII, we also consider the deposition of kinetic energy 
into galactic-scale winds, with initial velocities that vary linearly 
with the local velocity dispersion of the dark matter halo. The mass 
loading is then determined from the available energy from SNII ex¬ 
plosions (Springel & Hernquist 2003). With these choices for the 
sub-resolution physics, the Illustris simulations have been shown 
to reproduce several of the observational properties of galaxies, 
such as their stellar mass and color distribution, morphological va¬ 
riety, angular momentum distribution and typical scaling relations, 
among others (Torrey et al. 2014; Genel et al. 2014, 2015; Snyder 
etal. 2015; Sales et al. 2015). 

In the highest resolution run, or Illustris-1 (used here), the 
mass per DM and baryonic particle are moM ~ 6.26 x 10® Mq 
and mb ~ 1.26 x 10® Mq, respectively. Gravity is softened 
with a spline kernel (Hernquist & Katz 1989) and the softening 
length scales with redshift but is always comparable to or better 
than e = 700(1400) pc for the stellar (dark matter) particles. The 
most well-resolved gas elements can have extents as small as 48 
pc. 

Galaxies are identified as self-bound structures at each out¬ 
put using a combination of a Friends-of-Friends (FoF) algorithm 
(Davis et al. 1985) for groups and a subsequent application of the 
SUBFIND algorithm (Springel et al. 2001; Dolag et al. 2009). The 
centers are marked by the position of the particle with the minimum 
potential energy and galaxy properties are computed within twice 
the half-mass radius of the stars for each object. The time evolution 
of a given object is tracked by means of LHaloTree merger trees 
(Springel et al. 2005; Nelson et al. 2015) from a total of 136 output 
snapshots covering the redshift range 2 = 50 to 2 = 0. 
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Figure 1. Projected stellar maps for the 12 most massive galaxy clusters in lllustris (M 200 ^ 5 x 10^^ Mq). Individual “clumps” correspond to satellite 
galaxies inhabiting each of the clusters. Circles indicate the extent of the virial radius. 
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Figure 2. Left: Cumulative stellar-mass function of satellite galaxies for the simulated clusters in our sample. Individual clusters are shown with colored 
dashed lines and the median is indicated with solid black line/symbols. Our selection criteria for dwarf galaxies (M* = [3 X 10® — 1 X Mq) is 

highlighted by the gray area. We find a total of 1071 cluster dwarfs. Note the good agreement with the SDSS observed satellite mass function (taken from 
Yang et al. (2009), orange shading), particularly in the selection region of our dwarfs. Right: The color-stellar mass distribution of dwarfs in our cluster sample 
(red dots) and field control sample (blue crosses), selected to be central (not satellites) objects in the same mass range. Encouragingly, cluster dwarfs are much 
redder than the field sample, as indicated by observations. We show in grey circles/triangles SDSS results for galaxies in low/high density environments taken 
from Kauffmann et al. (2004). 


2.1 Simulated Dwarfs in Clusters and in the Field 

Our sample of dwarfs consists of a set of “cluster” dwarfs selected 
to inhabit the most massive groups and clusters in Illustris, and 
“field” dwarfs, a control sample consisting of central low mass 
galaxies that are isolated. The selection is done once at z = 0 
and the time evolution of both populations is studied by follow¬ 
ing the progenitors of each object through time. To construct the 
cluster sample, we consider all galaxy clusters with virial mass 
^^200 ^ 5 X 10^® Mq, where the virial radius r2oo is defined 
to be that which encloses 200 times the critical density of the Uni¬ 
verse, and the virial mass M 200 is the total mass contained within 
r 200 - Applying these criteria, we obtain a total of 12 clusters, with 
the most massive reaching M 200 = 2.32 x 10^'* M©, similar to 
mass estimates for the Fornax (Jordan et al. 2007a) and the Virgo 
(Bohringer et al. 1994) galaxy clusters. 

In Fig. 1 we show stellar density maps of the clusters in our 
sample, with the virial mass and original FoF-ID quoted in each 
panel. All clusters are populated by a large number of satellite 
galaxies, seen as individual “clumps” in each image, which cover a 
wide range of masses and properties. By definition, satellite galax¬ 
ies are all objects identified by SUBFIND that are not the centrals of 
their FoF group. 

The left panel of Fig. 2 shows the cumulative stellar-mass 
function of these satellite galaxies, indicated by colored dashed 
lines for each of the selected 12 clusters. We consider only satellites 
that are at 2 = 0 within the virial radius of their host system. The 
median of our sample (black symbols/solid line) agrees well with 
the satellite stellar-mass function of groups and clusters of similar 
mass in SDSS, as indicated by the orange shaded region (Yang et al. 
2009). Cluster dwarf galaxies are then selected from these satellites 


by a simple mass cut, 3 x 10® < M^/Mq < 1 x 10^° (gray re¬ 
gion), implying that all our dwarfs are resolved with at least 240 
stellar particles today. This selection criterion yields a total of 1071 
cluster dwarfs. 

For comparison, we also define a sample of “field” dwarfs in 
the same stellar mass range but that are central galaxies of their FoF 
group (i.e. they are not satellites of any larger system). A total of 
28391 dwarfs satisfy this condition, from which we randomly se¬ 
lect 1071 to match the number of cluster dwarfs. Encouragingly, 
this simple selection procedure results in interesting differences 
between the color distributions of the “cluster” and “field” dwarf 
populations, as expected from observations of galaxies in differ¬ 
ent environments. The right panel of Fig. 2 shows that at fixed 
stellar mass, the cluster dwarfs (red dots) are significantly red¬ 
der than the field objects (blue crosses), in good agreement with 
the observed redder colors of satellites around massive hosts (e.g.. 
Dressier 1980; Martinez et al. 2002; Wang et al. 2014). This sug¬ 
gests that our simulation properly captures the main physical pro¬ 
cesses driving the evolution of galaxies in different environments 
and can therefore shed light on the formation and evolution of dEs 
in cluster environments. 


3 THE ASSEMBLY OE DWARE GALAXIES IN 
DIFEERENT ENVIRONMENTS 


The cluster dwarfs sample a wide variety of orbits around their 
central hosts. We show several examples in the top row and first 4 
columns of Fig. 3. Dwarf satellites initially follow the Hubble ex- 
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Figure 3. Examples of the orbits (top row) and mass evolution (bottom row) of cluster dwarfs (Satellites A-D) and a field dwarf (rightmost column). Orbits 
show a wide diversity, with some dwarfs completing more than 3-4 revolutions around the cluster (example Sat. A) and others only recently arriving (Sat. 
D). The time evolution of the host cluster’s virial radius r 200 are indicated with a black dashed line in each panel. The time of maximum total mass (see 
green-dotted line in the bottom row) is indicated in all panels with a vertical dotted line. This time corresponds roughly to the moment they stopped being 
centrals to become satellites; which can happen right before crossing the virial radius of the cluster like in Sat. A, B and D, or before, if they were accreted into 
a smaller group first and then entered the cluster (like Sat. C). The mass evolution in the bottom row shows clear correlations with the orbits, with a decrease of 
mass after infall as well as close pericentric passages. Tidal stripping is not strong for the stellar component (solid red), although dwarfs in tightly bound orbits 
can experience significant stellar mass loss (example Sat. A). Small vertical arrows show the times at which each component (dark-matter, stars, gas) reaches 
its maximum. Stars continue to build up after infall, as shown by the shift between the dotted vertical line and the red aiTows. Note that for field dwarfs, which 
are not exposed to stripping, the mass in all components peaks only at the present time. 


pansion until they are captured by the gravitational field of a cluster 
and then turn around (seen as the moment of maximum distance 
from the cluster). They then begin the process of infall and settle 
into their orbits around the host potential. Some dwarfs in our sam¬ 
ple have fallen in more than 10 Gyr ago, completing several revo¬ 
lutions around the cluster’s center by the present (e.g. Satellite A). 
On the other hand, some objects have only recently been accreted 
into the cluster (Satellite D), crossing the host virial radius less than 
2 Gyr ago (the time evolution of the host virial radius is shown with 
a dashed line). From this perspective, the sample of cluster dwarfs 
is a complex group of objects whose properties are expected to be 
diverse in light of the variations present in their orbital histories. 

Orbits have a clear impact on the mass evolution of cluster 
dwarfs, as shown in the bottom row of Fig. 3. Dwarfs initially grow 
as expected (green dotted line shows the evolution of total mass) but 
they later reach a maximum after which their total mass decreases 
due to tidal stripping. The maximum is usually reached at the last 
time the dwarfs are the central object of their FoF group. Therefore, 
this time of maximum mass, shown as vertical dotted lines in each 
panel, is a good indicator of the time of infall, or the time when 
dwarfs stopped being central objects to become satellites. In what 
follows, we will refer to the time of infall finf or infall redshift Zint 
as the time when the simulated dwarfs reached their maximum total 

mass fniax- 

There are two interesting points to notice from these growth 


curves. Some dwarfs (e.g. Satellite C) first become a satellite of a 
different system that later fall into the cluster. This pre-processing 
of dwarfs is not uncommon, and is usually accompanied by a wig- 
gly pattern in their orbits as they move around a different host on 
the way to their final infall into the cluster (see Satellite C around 
t ~ 5.5 Gyr). Pre-processing effects in groups have been hypoth¬ 
esized to account for the observed properties of the Magellanic 
Clouds (Besla et al. 2007, 2010, 2012), and is also expected to hap¬ 
pen regularly in cluster dwarfs (e.g. Lisker et al. 2013). For our 
sample, 65% of the objects infall first into an intermediary system, 
where they spend a rather short amount of time (median of ~ 1.3 
Gyr), before finally infalling into their host galaxy clusters where 
they stay until 2 = 0. Second, the small difference between the 
green dotted and black solid lines suggests that stripping mainly 
affects the dark matter component (solid black line). The times 
when the maximum total mass fmax and maximum dark matter 
mass fmax are reached always track one another (see black arrows). 

Baryons, instead, show a different behaviour. Stars are much 
more centrally concentrated than the dark matter component, re¬ 
sulting in a small fraction of the stellar mass being tidally stripped 
(with some exceptions as in Satellite A). In our sample of 1071 
cluster dwarfs, 62% retain more than 80% of their maximum stel¬ 
lar mass and only 1% have lost more than 90% of their stars by 
2 = 0. For comparison, this corresponds to 2% and 18% objects 
respectively for the same thresholds in the dark matter component. 
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Figure 4. Distributions of the time of maximum total mass fmax (top) and 
maximum stellar mass (bottom). Histograms are normalized to the 

peak value in each sample (values quoted). Cluster/Field dwarfs are shown 
in red/blue. Due to stripping, dwarfs in clusters assemble their mass earlier, 
typically reaching their maximum mass at tmax ~ 5.5 Gyr or 2 : ~ 1, unlike 
the field population which continues to grow until today (tmax ~ 13.7 
Gyr). The time of maximum stellar mass in cluster dwarfs is about 3.5 Gyr 
later than the time of maximum total mass, indicating that satellite dwarfs 
can continue to form stars for a significant period of time after they have 
become satellites. 

Moreover, the stellar mass continues to grow after dark matter strip¬ 
ping has begun, as indicated by the red arrows marking the time of 
maximum stellar mass for each object, fmax- This agrees well with 
previous findings in Illustris of significant star formation continu¬ 
ing after satellite infall (Sales et al. 2015; Rodriguez-Gomez et al. 
2015). The gaseous component (dashed blue line) stops growing 
around the time of infall and tends to be equally or more greatly 
affected by the environment than the dark matter component. Tidal 
stripping as well as gas consumption from star formation and hy- 
drodynamical effects like ram-pressure stripping contribute to this 
effect. Most cluster dwarfs (74%) show significant gas depletion, 
losing more than 90% of their maximum gas mass, which is con¬ 
sistent with the low levels of star-formation driving the red colors 
seen previously in Fig. 2. 

The evolution of dwarfs in the field is remarkably different. 
The last column in Fig.3 shows the mass evolution of a random 
dwarf from the field sample. Of note, all components show a steady 
mass growth until 2 = 0, with occasional pronounced mass-gain 
events associated with mergers (for example at f ~ 4 Gyr). Galax¬ 
ies in the field do not experience significant stripping, and therefore, 
reach a maximum mass at the present day, as shown by the vertical 
dotted line and small arrows. This is not a feature of this particu¬ 
lar dwarf, but rather the norm in the whole field dwarf sample, as 
shown in Fig.4. The top panel shows the distribution of fmax (the 
times when the total maximum mass is reached) for all cluster (red) 
and field (blue) dwarfs. The median and 25%-75% percentiles are 
shown as vertical dashed and dotted lines, respectively. For most 
dwarfs In the field, fmax ~ 13.7 Gyr (i.e., ‘today’) whereas the 
median for cluster dwarfs is tmax ~ 5.6 Gyr or 2 = 1.1. A sim¬ 
ilar behaviour is seen for the time of maximum stellar mass, with 


redshift 

6 3 2 1.3 1 0.5 0 



Figure 5. Median stellar mass evolution for the sample of cluster (red) and 
field (blue) dwarfs. The 25%-75% dispersion on the median is shown as the 
shaded area for the cluster dwarf sample only, but is of similar magnitude 
for the field dwarfs. Although both samples have (by construction) the same 
mass at 2 = 0 their time evolution is significantly different: cluster dwarfs 
were more massive than the field sample at earlier redshifts. Both samples 
gain mass at approximately the same rate up to 2 ~ 1, when the growth 
of cluster dwarfs decelerates due to environmental effects (median infall 
time of the sample shown by the vertical dashed line). In contrast, a sample 
of dwarf analogs that had the same properties as the progenitors of cluster 
dwarfs at the time of infall (instead of at z=0) would outgrow the median 
stellar mass in cluster dwarfs by up to a factor ~ 3 today (green line). 

field dwarfs assembling their stars significantly later than their clus¬ 
ter counterparts (f^ax ~ 13.6 and 9.14 Gyr for field and cluster 
dwarfs respectively). Fig. 4 highlights that a significant level of star 
formation after infall is a common feature for cluster dwarfs (the 
median fmax and f^ax differ for more than ~ 3 Gyr); a subject we 
will return to in Sec. 4. 


This different stellar mass assembly can be better appreci¬ 
ated from Fig. 5 for the population of cluster (red) and field (blue) 
dwarfs. Curves correspond to the median of the samples, and the 
25%-75% dispersion is indicated by the shaded region (for clarity 
only shown for the cluster dwarfs but is of similar magnitude for the 
field objects). Because of our selection, the curves overlap at 2 = 0 
where both samples were chosen to be within the same stellar mass 
range: 3 x 10® < M^/Mq < 1 x 10^°. Flowever, going back 
in time, the samples diverge. Field dwarfs form their stars steadily 
throughout their history up to 2 = 0. In contrast, dwarfs in clusters 
grow their mass rapidly at early times, then being 3-4 times more 
massive than the field sample. However, by the median time of in¬ 
fall of the cluster sample, (finf) ~ 5.5 Gyr (vertical dashed line), 
the stellar mass freezes and does not appreciably change thereafter. 
This behaviour is caused by more dwarfs becoming satellites, low¬ 
ering their star-formation rates and being exposed to tidal stripping. 

This excess of mass at early times with respect to today’s field 
objects suggests that, if they were in the field, cluster dwarfs would 
have grown significantly more in mass by today. Processes asso¬ 
ciated with the cluster environment have truncated that growth. 
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Figure 6. Median star formation rate (SFR, top panel) and the specific star 
formation rate (sSFR, bottom panel) as a function of time for cluster (red) 
and field (blue) dwarfs. Color coding is the same as in Fig. 5. Dwarfs in 
clusters form their stars at higher SFR than the field sample. Flowever, the 
difference is only due to their larger stellar mass at high redshifts, since 
the sSFR of both samples agrees well before the average time of infall 
(finf) ~ 5.5 Gyr. After that, the median rate of star formation in cluster 
dwarfs decreases, with most objects becoming quiescent by z = 0. We de¬ 
fine sSFR = 10“^^ as our threshold to divide star-forming from 

quiescent objects (black dashed line). 


We examine this possibility by defining a group of cluster dwarf 
“analogs” (green line), selected to have similar mass as the cluster 
dwarfs at time of infall but that, unlike the cluster dwarfs, remain 
isolated field objects until the present^. This exercise shows that 
a sample selected to match cluster dwarfs at infall rather than at 
2 = 0 ends up a factor ~ 2.5 — 3 more massive on average than 
the cluster dwarfs today. The most appropriate comparison sample 
is thus unclear and may depend on the question being addressed. 
Our results indicate that selecting isolated objects in the field hav¬ 
ing ~ 3 times more mass than the cluster dwarfs we are targeting 
may provide a better match to the early phases of growth of both 
samples. To facilitate the comparison with observations, in what 
follows we will still consider our “field” objects selected within 
the same mass range as cluster dwarfs at 2 = 0 to be the control 
sample, but we bear in mind the lesson learned from the “analogs” 
sample. 


^ In practice, the selection of “analogs” is done as follows: at the median 
infall time of the cluster dwarfs sample (z) ~ 1.1, we select all field objects 
that at that time have a stellar' mass within a factor of 2 of the median stellar 
mass of the progenitors cluster dwarfs at the same time. We then follow 
the selected objects forward in time to identify those that remain central 
galaxies to z = 0 and choose 1071 objects randomly from that set. 


4 THE STAR FORMATION HISTORY OF DWARF 
GAFAXIES: CFUSTER VS. FIFED 

We now turn our attention to how dwarfs form their stars according 
to their environment. Fig. 6 shows the median star formation rate 
(SFR, upper panel) and specific star formation rate (sSFR, lower 
panel) of our samples as a function of time. We define the sSFR 
as the SFR divided by instantaneous stellar mass of the galaxy. As 
before, red/blue corresponds to cluster/field dwarfs and 25%-75% 
percentiles are indicated by the shadowed region. As expected from 
the mass difference seen in Fig. 5, cluster dwarfs form stars at a 
higher rate than in the field at early times, but the trend reverses as 
more dwarfs become satellites after 2 = 1. By the present day, the 
vast majority of cluster dwarfs have stopped forming stars, consis¬ 
tent with the red colors displayed in Fig. 2. Interestingly, the SFR 
difference at high redshifts is only due to the different masses of the 
samples: once the SFR of each dwarf is normalized by their (instan¬ 
taneous) M*, the specific rates between the samples are consistent 
with each other (bottom panel). Cluster dwarfs are not “special” ob¬ 
jects in this regard and simply form stars according to the expected 
rate given their stellar mass. 


Beyond averaging the samples, examination of individual star 
formation histories of dwarfs reveals another interesting aspect of 
cluster dwarfs. Fig. 7 shows the SFR evolution corresponding to 
Satellites A-D and the field dwarfs shown in Fig. 3. The curves 
in the first three panels show significant peaks or “bursts” of star 
formation that are not present in the field objects or in Satellite D. 
Such bursts, when present, seem to occur around the time of infall 
(indicated by the vertical dotted line) and we have checked that 
these events are usually associated with either the crossing of the 
virial radius of the cluster or to the time of first pericentric passage. 
Often times the star formation is quenched immediately following 
this final “burst”. 

We would expect that the tidal and ram-pressure interac¬ 
tions with a cluster could trigger a compression of the interstellar 
medium of a cluster dwarf, temporarily raising its star formation 
rate. Such an effect has previously been seen in idealized simula¬ 
tions (Barnes & Hemquist 1996; Mihos & Hernquist 1996; Mayer 
et al. 2001). Here we find that the triggering of short “bursts” of 
star formation by the environment should be quite common among 
dwarfs in clusters today. Additionally, we have explicitly checked 
that mergers with other galaxies within the clusters are far too in¬ 
frequent to be responsible for the measured star formation increase 
in cluster dwarfs. At any given time, satellite dwarfs have a median 
merger rate at least ~ 1 dex below that of the field counterparts. 

We quantify the presence of bursts of star formation by mea¬ 
suring the standard deviation of the SFR at each snapshot with re¬ 
spect to the time-average: 


- SFR)VNt 


where Nt is the number of time points where we measure the SFR 
(i.e., the number of output snapshots of the simulation), and SFR 
is the time averaged SFR from all snapshots (indicated with a cyan 
curve in each panel of Fig. 7). Values for ctsfr are quoted for all 
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Satellite A Satellite B Satellite C Satellite D Field 



Figure 7. Examples of individual star formation histories for the same satellite and field dwarfs in Fig. 3. As in the case of the orbits, the SFR histories of 
cluster dwarfs exhibit large variations. It is not uncommon to see “peaks” or episodes of intense star formation (for example Sat. A, B) that coincide with the 
crossing of the virial radius of the cluster and/or first pericentric passages. Starburst events are not common for field dwarfs. Some cluster dwarfs show a more 
constant SFR that is comparable to field objects (for example Sat. D). In each panel, we quote the “starburstiness” of the curve, itsfr, defined as the standard 
deviation of the SFR history with respect to the time-average (crspj^) (cyan line). Large values of (cspr) indicate the presence of starbursts. (See text for 
more details.) 


dwarfs in Fig. 7 and they appear to correlate well with the visual 
intuition of “hurstiness” in each curve^. 

Fig. 8 extends the trend hinted hy these 5 individual exam¬ 
ples to all dwarfs in our samples. We show ctsfr for cluster (red & 
orange dots) and field (blue crosses) dwarfs, plotted as a function 
of the local density 5 (where S is defined as the galaxy overden¬ 
sity within the radius that encloses the 5th nearest neighbour, and 
only galaxies with r-band magnitude brighter than —19.5 are used 
to compute the density field). Not only do field and cluster dwarfs 
separate cleanly in environment (as expected), they also show very 
different distributions of ctsfr- Cluster dwarfs typically have larger 
CTSFR values, indicating the presence of significant peaks or star- 
burst events in their star-formation histories. 

Moreover, green dots/histogram in Fig. 8 show the “analogs” 
sample described in Sec. 3 (selected to have similar progenitor 
masses as cluster dwarfs at time of infall). This sample is also char¬ 
acterized by low CTSFR values alike the field galaxies, confirming 
that SFR enhancements in cluster dwarfs are purely associated to 
environmental effects and is not a bias in their progenitor mass. Al¬ 
though the trend with environment is weak above the field-cluster 
threshold, dwarfs in the very inner regions of clusters show even 
larger ctsfr compared to the rest of the satellites (the orange his¬ 
togram includes only satellite dwarfs that are today within 0.5r2oo)- 
Our results indicate that cluster dwarfs, particularly those in the 
inner parts of clusters today, likely had a more “bursty” star for¬ 
mation history than dwarfs in the field. 

Star formation in satellite galaxies is expected to be trun¬ 
cated sometime after infall. Because the time-scale for this to oc¬ 
cur sheds valuable information on the mechanisms responsible for 
this quenching (and how it may depend on host and satellite prop¬ 
erties), the subject has received significant attention from the com¬ 
munity (e.g., Wang et al. 2007; Font et al. 2008; Wetzel et al. 2013; 

^ We have explicitly checked the robustness of itsfr by first smoothing 
the curves with a Savitzky-Golay filter of width t = 0.5 — 1 Gyr. Since 
all measures seem to correlate well with the unfiltered version we decided 
to use the data without the smoothing, since the filtering introduces an arbi¬ 
trary time-scale that is in principle unknown and could vary from dwarf to 
dwarf. 


Wheeler et al. 2014; Fillingham et al. 2015). In Fig. 9 we consider 
the time elapsed between infall and “quenching”, defined as the 
first time a dwarf’s SFR drops below sSFR ^ 10“^^yr“^. This 
threshold was chosen to be compatible with the definition in Wet¬ 
zel et al. 2013, and we have additionally confirmed that it applies 
to the scale of dwarfs by checking that it is always 2 dex below 
the average sSFR of dwarf centrals at any given time. For cluster 
dwarfs, we show the quenching time as a function of satellite mass, 
where dots correspond to our simulated dwarfs that are quiescent 
by 2 = 0 (cluster dwarfs that are still star-forming today do not 
appear in this plot; they amount to ~ a third of the sample). 

Interestingly, we find a positive trend with mass: the quench¬ 
ing timescale increases with the stellar mass from a median of ~ 3 
Gyr for our smallest objects to ~ 5.5 Gyr at M* = 10^° Mq (solid 
dark red curve). Notice that the scatter about this mean trend is sig¬ 
nificant, comfortably spanning the range ~ 1-8 Gyr after infall. 
Moreover, the dispersion is not random, but appears to correlate 
strongly with the mode of star formation, as quantified by ctsfr: 
dwarfs with significant “hurstiness” (large ctsfr values) prefer¬ 
entially show shorter time-scales for quenching than objects with 
a lower ctsfr at a given stellar mass (see color-coding). Long 
quenching time-scales are associated with dwarfs that passively 
consume their gas while rapid quenching is often triggered by a 
last burst of star formation that consumes all available fuel for the 
formation of stars. 

Comparing these results to observational data is not straight¬ 
forward. Samples of groups/clusters are typically not complete 
down to the masses of the dwarf galaxies studied here. On the other 
hand, the properties of dwarf galaxies in this mass range are well- 
studied in the Local Volume, where they are easier to observe due 
to their faint luminosities, or around ~ L, isolated primaries where 
statistics are better for the hosts compared to more massive groups 
and clusters. Although not ideal, we attempt to put our results in 
the context of available observations in Fig. 9. 

We start by comparing with the quenching time-scales in clus¬ 
ters but for more massive galaxies than our dwarfs (where there is 
observational data). For this, we extended our analysis to all satel¬ 
lite galaxies in the 12 simulated cluster hosts in Illustris. We select 
all surviving satellites at 2 = 0 that are within the virial radii of our 
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clusters and define their quenching time-scales in the same way as 
earlier. The median values obtained are indicated with the dark red 
dashed curves in Fig. 9. For such massive satellites, we overplot 
observational estimates from Wetzel et al. (2013) using SDSS data. 
The agreement is good, as shown by the solid red/green curves that 
bracket our host halo masses. Unfortunately, this does not guaran¬ 
tee that the time-scales for lower mass dwarfs in the simulations 
is reasonable; but it is reassuring that the model is similar to the 
observational data in the mass range where they overlap. 

A closer look at Fig. 9 reveals a preferred mass around M* ~ 
8 X 10® Mq where this quenching timescale is maximum, indicated 
by the medians in solid and dashed dark red curves. The reason for 
this upturn is unclear, but is likely connected to the physical mecha¬ 
nisms quenching these galaxies; we defer to a detailed study of this 
in a companion paper (Mistani et al., in-prep). The timescales for 
our M* ~ 10® Mq objects seem short compared to similar mass 
objects observed around L* galaxies (data taken from SDSS data 
in Fillingham et al. (2015) and Wheeler et al. (2014)); however the 
comparison is inconclusive, as we cannot rule out a dependence on 
host halo mass on this scales. A more detailed comparison awaits 
until observational constraints on the quenching of dwarfs in clus¬ 
ters become available. 

Summarizing, our results indicate that the star formation 
histories of cluster and field galaxies -with similar stellar mass 
today- differ; cluster dwarfs form their stars earlier, at higher star- 
formation rates and with occasional starbursts as they orbit their 
host cluster. In contrast, field dwarfs form their stars at later times 
and with typically lower SFRs. This clearly leaves an imprint in the 
stellar populations of these two sets of dwarfs, explaining the much 
redder colors in the cluster dwarfs compared to objects in the field. 
Could these different star formation histories also leave an imprint 
in the globular clusters populations of dwarfs? We investigate this 
in the following Section. 



1 + 5 Number 


Figure 8. Measurement of “starburstiness” in SFRs of dwarf galaxies 
(quantified by means of (tsfr) as a function of the galaxy density of the 
local environment 5. Large (tsfr are associated with the presence of star- 
bursts. Field dwarfs (blue crosses) inhabit less dense environments and 
show lower ctsfr values compared to cluster dwarfs (red circles). The early 
“analogs” of cluster dwarfs (green symbols) also show low (tsfr consistent 
with field objects. Flistograms on the right show the distribution of crgpjj^ for 
both samples together with the median indicated by dotted lines. Starbursts 
are more common among cluster dwarfs, particularly those that are today 
in the densest inner regions of clusters (orange distribution corresponds to 
cluster dwarfs within r < 0.5r2oo today). 


5 IMPLICATIONS FOR THE EXCESS OF GLOBULAR 
CLUSTERS IN DWARE ELLIPTICAL GALAXIES 

There is not a clear picture for how globular clusters (GCs) form, 
but their properties are consistent with being the evolved versions 
of the young star clusters that are found in nearby star-forming 
galaxies (Ho & Filippenko 1996; Mengel et al. 2002; de Grijs et al. 
2004). For example, the GC population of dwarfs in the Local Vol¬ 
ume appears consistent with the initial mass function of young stel¬ 
lar clusters, after accounting for stellar and tidal evolution (Krui- 
jssen & Cooper 2012). Therefore we consider the possibility that 
the processes and physical conditions behind the formation of to¬ 
day’s GCs are the same than those driving the formation of young 
star clusters at 2 : = 0. 


Analytical arguments suggest that gas-rich, high-pressure en¬ 
vironments offer the most propitious conditions for the forma¬ 
tion of bound star clusters (Harris & Pudritz 1994; Elmegreen & 
Efremov 1997; Shapiro et al. 2010; Kruijssen 2012); a claim that 
has also found support from numerical simulations (Kravtsov & 
Gnedin 2005; Bournaud et al. 2008; Renaud et al. 2008). Observa- 
tionally, there is strong evidence in favour of a link between star 
cluster formation and local properties of the interstellar medium. 
In particular, the fraction of stars born in clusters (that survive the 


embedded phase) compared to the total amount of stars formed, 
termed F, appears to correlate with the density of star formation as 
(Goddard et al. 2010); 


r 


0.29 


SsFR 

Mqui — ^kpc~^ 


0.24 


( 2 ) 


with Esfr defined as the averaged projected density of star forma¬ 
tion in a galaxy. Such a scaling was later confirmed by other in¬ 
dependent studies (Adamo et al. 2010; Silva-Villa & Larsen 2011; 
Adamo et al. 2011). If indeed star clusters form at the high den¬ 
sity tail of the interstellar medium (Kruijssen 2012), then it is ex¬ 
pected that GC formation was more efficient at higher redshifts 
(when galaxies were typically more gas-rich and turbulent) and 
during starburst events; interestingly this corresponds to two of the 
main differences in the star formation histories of our cluster dwarfs 
compared to our field galaxies. 

To address this issue in some detail and to roughly predict the 
GC efficiencies in cluster and field dwarfs, we implement a simple 
model for the formation of GCs as a post-processing algorithm in 
our simulations. The model involves three steps; i) compute the 
amount of stellar mass born in star clusters, ii) sample an initial 
mass function of clusters and Hi) use a particle tagging scheme to 
address the tidal removal of GCs from dwarfs orbiting inside galaxy 
clusters. 

For the first step, we use the empirical relation from Goddard 
et al. (2010), where for each dwarf at a given time we compute 
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Figure 9. Quenching time-scales for dwarfs in clusters as a function of 
stellar mass, where the time-scale is defined as the difference between in¬ 
fall time and the time at which the sSFR drops below 10~^^ yr“^. Clus¬ 
ter dwarfs are shown with circles color coded by their crsFR- The median 
quenching time increases with stellar mass (red solid line). Interestingly, 
the dispersion around this median is large but is well-correlated with the 
mode of stai* formation: starbursts events (large ctsfr) associated pref¬ 
erentially with the shortest timescales at a given stellar mass. To compare 
with observations, we extend this analysis to satellite galaxies in clusters 
but that are more massive than our dwarfs (median and dispersion shown 
with red dashed lines). These agree well the time-scales derived for SDSS 
galaxies in Wetzel et al. 2013 (red/green lines) for this mass range. The 
black square shows the inferred quenching time-scale for 10^ Mq dwarfs 
orbiting around L* isolated primaries. 
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Figure 10. Mass in GCs predicted for the field (blue) and cluster (red) 
dwarf samples according to our post-processing model of GC formation 
(see Sec. 5). Although the distributions overlap, dwarfs in clusters have on 
average ~ 2 times more mass in GCs than the field dwarfs; an effect domi¬ 
nated by the earlier and more bursty star-formation in cluster dwarfs. 


In order to compare more closely with observational measure¬ 
ments of the GC efficiency in dwarfs, we need to translate our es¬ 
timates of mass in GCs to GC numbers. Following Jordan et al. 
(2007b), we assume a Schechter-like cluster mass function and ran¬ 
domly sample from such a distribution until we account for the total 
GC mass for each dwarf. In particular, we adopt: 


the fraction of stars formed that are born in bound clusters using 
Eq. 2. We define Esfr as the sum of the star formation of all gas 
cells in the galaxy and we divide by the area of a circle with radius 
equal to that enclosing 95% of the star-forming gas (a good proxy 
for the Ha radius). We then repeat this calculation throughout the 
time evolution of each dwarf, averaging in time-steps of 1 Gyr (we 
have explicitly checked that the results do not depend on the time 
window used). Adding up the mass in stars that is bom in clusters 
for each of the time-steps gives us a prediction for the total mass 
expected in star clusters for each dwarf. 


Fig. 10 shows the distribution of the computed mass in GCs 
for the population of cluster (red) and field (blue) dwarfs. Although 
the distributions overlap, we find a significant difference between 
the two, with cluster dwarfs having on average a factor ~ 2 times 
more stellar mass in bound clusters than field dwarfs (vertical lines 
show the median in each sample). The difference is largely due to 
the earlier and higher SFRs in the former, as discussed in Sec. 4. 
This overall trend is in good agreement with results from the semi- 
analytical models presented in Peng et al. (2008), but the star- 
formation histories of field and cluster dwarfs are much more alike 
in their model (based on the Millennium simulations plus Guo et al. 
(2011) semi-analytical catalog) than in our sample. 


— cxM-^exp{-M/M,), (3) 

where Me = 2 x 10® Mq is the truncation mass. For our study, 
we focus on the high mass end, where GCs are more resilient to 
destmetion by internal and environmental mechanisms. We sample 
globular clusters in a mass range of [1 x 10*Mq, 1 x 10 ®Mq], 
where the upper bound is set to ensure that GCs have less mass 
than the galaxy itself. This is in general not important, as the rapid 
fall-off of the Schechter function ensures that high-mass clusters 
are highly uncommon. Notice that the low mass end we choose for 
the sampling works as a “normalization” factor to our total num¬ 
ber of GCs (if we choose a smaller mass, we get a higher number 
of GCs and the same mass can now be distributed among many 
smaller-mass GCs than before). Thus, in what follows, our results 
are differential between both the dwarf samples. Our main aim is 
to compare the results for the cluster and field samples with respect 
to each other under the same assumptions. 

For dwarfs in the field, this is enough to compute their GC 
specific frequency Tn, defined as the number of GCs divided by 
the stellar mass of the dwarf, Tn = Ngc/Mt/lO^ Mq (Zepf & 
Ashman 1993). However, for dwarfs in clusters, the tidal forces ex¬ 
erted by the host cluster potential can potentially strip some of the 
GCs originally formed in-situ. We therefore implement a particle¬ 
tagging technique, where we “paint” GCs on some dark-matter par¬ 
ticles and use those particles as tracers of the time evolution for the 
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Figure 11. Final specific frequency of GCs for the sample of dwarfs in clus¬ 
ters (colored circles) and in the field (black crosses). Tjv is computed as the 
number of GCs within an aperture of 8 kpc divided by the stellar mass of 
each object. Tjv also accounts for the effects of environmental tidal strip¬ 
ping. In agreement with observations, we find that cluster dwarfs have a 
larger frequency of GCs than field dwarfs of the same , a trend that is 
enhanced for dwarfs with stai'bursty star formation histories (higher crspjc, 
see color coding) and/or dwarfs in the inner regions of clusters. Magenta 
and cyan lines show the medians for cluster dwarfs with r < 0.3r2oo and 
r > 0.7r200, respectively; field dwarfs shown with dotted black line. (For 
clarity, we only show 1 in 5 field dwarfs selected randomly from the sam¬ 
ple.) Inclined grey lines show lines of constant number of GCs: 100,10 and 
1 from top to bottom. Differences are relative between the field and cluster 
samples, as the vertical normalization is chosen to coincide roughly with 
observed Tjv range (see text for more detail). 


GC population. This technique is not dynamically self-consistent 
(the mass of the dark-matter particle does not coincide with the GC 
that was “painted” onto it) but the expected error is small. The tech¬ 
nique has been successfully used before in studies of GCs (Griffen 
et al. 2010; Ramos et al. 2015), stellar halos (Cooper et al. 2010, 
2013), brightest cluster galaxies (Napolitano et al. 2003; Laporte 
et al. 2012; Cooper et al. 2015), and structure of dwarf spheroidal 
galaxies (Penarrubia et al. 2008), among others. 

For our tagging technique, we randomly select particles fol¬ 
lowing a steep density profile, poc oc which is consistent 

with the distribution of GCs around the Milky Way according to the 
GC catalog in Flarris (1996). We use a maximum distance equal to 
10 kpc, but the steep slope ensures that the half-number radius is 
only 1-2 kpc, in agreement with observations of dwarfs (Georgiev 
et al. 2009). Our results do not depend on this maximum radius 
assumed. For each cluster dwarf, the tagging in time is done only 
once, at fmax- Although some of the star clusters will continue to 
form after infall into the galaxy cluster, we assume that the total 
number of GCs is already in place at the time of infall, which sim¬ 
plifies the procedure compared to tagging particles at individual 
time-steps. This approach is conservative, as tagging particles at an 
earlier time can only increase the chances of stripping later on. 

After the particles have been selected at fmax for each dwarf, 
we use their IDs to check for their positions at 2 = 0 and consider 
the GCs as bound/stripped according to their distance to the dwarf 


where it was bom. In order to approximate the aperture used in 
studies of the Virgo cluster (Jordan et al. 2007b, Peng et al. 2008), 
we use Tcut = 8 kpc, but have also checked that apertures in the 
range 5-10 kpc yield similar results. For field dwarfs (not exposed 
to stripping) the tagging is not necessary, and we simply do the spa¬ 
tial sampling once at fmax and check for the number of GCs within 
the aperture. In order to average-out the results of a particular ran¬ 
dom realization, we repeated the tagging exercise 10 times for each 
dwarf and assigned the final Tjv based on the average of these 10 
realizations. 

Fig. 11 shows the final specific frequency of GCs, Tjv, pre¬ 
dicted by our model using all massive GCs with mass above M — 
4 X 10® Mq (this mass cut-off is not important, but was chosen 
to agree roughly with the average normalization of observed GC 
frequencies in dwarfs). We show the predicted T^v as a function 
of stellar mass for the sample of cluster (solid colored dots) and 
field (black crosses) dwarfs. We find a clear excess of GCs for clus¬ 
ter dwarfs, even after accounting for tidal stripping within the host 
clusters, with the median Tjv of cluster dwarfs about ~ 1.65 times 
larger than dwarfs in the field. Most importantly, we find a clear 
correlation between this GC excess and the “starburstiness” of the 
star-formation histories in the dwarfs, as shown by the color coding 
of the dots: cluster dwarfs with clear peaks in their star formation 
histories (larger ctsfr) preferentially dominate the tail of high GC 
frequencies, whereas dwarfs with a more quiet star formation his¬ 
tory tend to show Tm closer to those of field dwarfs. We therefore 
see that the combination of early star formation plus the starburst 
events associated with infall and pericentric passages for cluster 
dwarfs can account for the more numerous GC populations ob¬ 
served for dEs. 


Observationally, much of the excess in GC specific frequency 
comes from nucleated dEs whereas non-nucleated dEs have GC 
frequencies more similar to those of dirrs (Miller et al. 1998; Miller 
& Lotz 2007). Although our numerical resolution is not sufficient 
to address the morphology of the simulated dwarfs, the correlation 
between high Tjv and the presence of starbursts suggests that nu¬ 
cleated dEs are simply the tail of the distribution of dwarf galaxies 
where environment caused the strongest bursts of star-formation 
associated with infall and pericentric passages, in agreement with 
the proposal of Miller & Lotz (2007). Moreover, the excess in Tjv 
found in Eig. 11 is dominated by dwarfs located in the inner re¬ 
gions of the clusters. The magenta line shows the median specific 
frequency of GCs for dwarfs with r < 0.3r2oo which is a factor 
~ 1.5 larger than that of cluster dwarfs in the outskirts r > 0.7r2oo 
(cyan curve) and dwarfs in the field (dashed black curve). This ra¬ 
dial dependence is in agreement with the observed trend for dEs in 
the Virgo cluster (Peng et al. 2008). 

The effects of tidal stripping can be significant, but is not 
enough to erase the initial difference (Eig. 10) expected from our 
calculation of the mass in GCs. Eig. 12 shows a more comprehen¬ 
sive view of the final distribution of GCs (red dots) with respect to 
the dwarfs (green squares) within a galaxy cluster (in this case, the 
most massive of our host clusters). The image shows a luminosity- 
weighted stellar mass map and indicates that although the majority 
of GCs are still bound to the dwarfs in which they were bom, there 
are several that have been stripped by the cluster potential and may 
now belong to the background population of GCs of the massive 
central galaxy in the cluster. Dwarf elliptical galaxies seem then 
able to contribute a fair number of “accreted” GCs to the central 
potential of the host cluster. Moreover, tidally disrapted and unre- 
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Figure 12. Final distribution of tagged GCs (red dots) overlaid on the projected stellar-mass map of our most massive cluster. Dwarf galaxies in this cluster 
(the only contributors of GCs in this figure) are highlighted with green squares. Although most of the GCs remain bound to the dwarfs where they formed, 
several of them have been tidally stripped from their native dwarf galaxy, particularly in the inner regions of the cluster where tidal forces ai'e stronger. The 
stellar content and phase-space properties (positions and velocities) of such stripped GCs could still retain information of their dwarf progenitors, such as 
orbits, star-formation histories and also help determine the tidal radii of cluster dwarfs. 


solved GCs will contribute to the diffuse intra-cluster light instead, 
a study of which we defer to future work. 

The results of this section are encouraging, but we emphasize 
the limitations of our model. One strong assumption we have made 
is that the star cluster efficiency F scales, at all redshifts and for 
low mass galaxies as the one considered here, with the density of 
star formation as indicated by Eq. 2. The equation is a power-law fit 
to observations in nearby star-forming galaxies, and it is not guar¬ 
anteed that the same relation holds at higher redshifts. In order to 
explore alternative parametrizations of F, we considered the theo¬ 
retical models presented in Kruijssen (2012), which are based on 
physical properties of the gas and therefore can be applied at any 
given time. In this model, F scales with the surface density of the 
gas Eg, which we can measure from our simulations. 

Although the Kruijssen model is more detailed than our nu¬ 
merical resolution allows (it requires measurements for the stability 


parameter Q for galaxy disks, among other factors), using standard 
values of Q = 1.5, we were able to reproduce the factor ~ 2 rela¬ 
tive excess in GC mass for dwarfs in clusters compared to the field^. 
Additionally, we have checked that our definition of size used to 
compute Esfr in Eq. 2 does not change the results we presented. 
To do this we vary the definition of size using: the radius enclosing 
95% of the star forming gas (the proposed model), half-mass radius 
of the stars and half-mass radius of the gas, and in all cases recover 
the relative difference between cluster and field dwarfs presented 
in Fig. 11. 

A second important assumption in our model is that the de¬ 
struction of massive GCs is similar for dwarfs in clusters and 


® We thank the author for making the code to compute GC efficiency pub¬ 
licly available at http://www.mpa-garching.mpg.de/cfe. 
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in the field. However, for massive GCs, dynamical friction (DF) 
forces can accelerate their orbital decay, eventually causing them 
to merge with the host galaxy. Analytical estimates of time needed 
for GCs to sink to the center of their hosts due to DF can vary 
from 1-10 Gyr depending on GC mass for a host dwarf galaxy with 
M* ~ 10® Mq (see fig. 19 in Turner et al. 2012). This indi¬ 
cates that the magnitude of the GC excess in cluster dwarfs can 
be overestimated in our model. Nonetheless, the net effect of DF 
is far from clear. For example, the DF timescale depends strongly 
on the (unknown) initial distribution of the GCs and also on the 
dark matter density profile of the dwarfs. Halos with lower cen¬ 
tral densities or constant dark matter cores show much longer DF 
timescales (Hernandez & Gilmore 1998; Cole et al. 2012). Because 
tidal stripping of the dark matter halo in satellite dwarfs will tend 
to lower the central densities (e.g. Hayashi et al. 2003), the sink¬ 
ing of GCs in cluster dwarfs could result, overall, less efficient than 
in the field. Therefore, although some of the most massive GCs in 
cluster dwarfs may merge with their hosts, due to the steep mass 
function of GCs favoring in number lower mass GCs, the GCs fre¬ 
quency in cluster dwarfs is likely to remain significantly enhanced 
with respect to the field. 

Finally, we also tested the effects of assuming a shallower GC 
radial distribution (and thus more prone to tidal stripping for cluster 
dwarfs). Following observational results by Beasley et al. (2009), 
we assumed p oc (instead of the steeper —3.5 slope found in 
the Milky Way) and sample the position of GCs within a maximum 
radius of 5 kpc. This cut-off was chosen such that the radius con¬ 
taining half of the GCs is ~ 1-2 kpc as observed for dwarfs. We 
repeated the particle tagging step and check for the number of GCs 
that remained bound to their host dwarfs, as before. We find that al¬ 
though some cluster dwarfs are stripped more easily of their GCs, 
the majority of the sample manages to retain most of their GCs, 
suggesting that our result in Fig. 11 is robust to changes in the ra¬ 
dial distribution of GCs that are within observational constraints 
(see also Smith et al. 2013). 

We conclude that the model is not overly sensitive to the arbi¬ 
trary choices we have adopted and that it represents an attractive 
tool for studying the formation of GCs in cosmological simula¬ 
tions of galaxies where the resolution is not adequate to describe 
the scale of individual star clusters. The application of this model 
to our sample of dwarf galaxies indicates that dwarf ellipticals in 
clusters would naturally have an enhanced frequency of GCs com¬ 
pared to field dwarfs, an effect that is most prominent for dwarfs in 
the inner parts of clusters; and that such a difference is not only a 
reflection of an early assembly but is also due to the starburst events 
associated with infall into the host cluster. 


6 SUMMARY & CONCLUSION 

We use the Illustris simulation to study the formation of dwarf 
galaxies in clusters. We select 12 host galaxy clusters in the virial 
mass range M 200 = 0.72-2.32 x 10^^ Mq, with the high mass end 
akin to systems like the Virgo and Fornax clusters. Dwarf galaxies 
are selected to be within the virial radii of these clusters at 2 = 0 
and have stellar masses in the range M, = 3 x 10®-1 x 10^® Mq. 
We find 1071 candidate dwarfs that satisfy our selection criteria. 
We construct an equal-number control sample of field dwarfs se¬ 
lected in the same stellar mass range at 2 = 0 but that are central 
objects of their FoF group. Our main findings can be summarized 
as follows: 

• At 2 = 0 the population of dwarfs in clusters has redder colors 


and lower star formation rates than field objects, in good agreement 
with the properties of observed dEs. The parameters in the simula¬ 
tion were not tuned to reproduce this, so this is a successful predic¬ 
tion of our model that motivates further studies into the evolution 
that shaped these two different populations. 

• Dwarfs in clusters are mostly genuine low mass objects and 
not the descendants of more massive galaxies that have been sub¬ 
stantially stripped. Only 14% of our dwarfs have lost more than 
50% of their stellar mass due to tidal effects. 

• Although by construction both samples have the same stellar 
mass at 2 = 0, the assembly of cluster and field dwarfs is remark¬ 
ably different: cluster dwarfs form their stars at earlier times and 
at higher star formation rates compared to field dwarfs. Because of 
this difference in the build-up times, the progenitors of dwarfs that 
are in clusters today were, at the average infall time of the sam¬ 
ple 2 ~ 1, about 3.5 times more massive than the progenitors of 
today’s field dwarfs at the same time. 

• The star formation histories of dwarfs in clusters often show 
significant starburst events which are not present in the population 
of field dwarfs. These starbursts are associated with infall and/or 
first pericentric passages, likely due to gas compression in the 
dwarfs under the combined effects of tidal and ram-pressure forces. 

• The cessation of star formation in cluster dwarfs occurs within 
timescales (median) ~ 3-5.5 Gyr after they become satellite galax¬ 
ies. We find a positive trend with mass, in which low mass dwarfs 
go quiescent most rapidly. Nevertheless, there is a large scatter 
around these median values, with differences ^ 7 Gyr at the 
same M*. This highlights the varied paths by which star formation 
is halted: dwarfs with significant starburst events quench rapidly, 
whereas others in the sample simply consume their gas fuel and 
quench over much longer timescales. 

• The different star formation histories of cluster and field 
dwarfs can have a significant impact on their populations of globu¬ 
lar clusters. To study this, we implement a post-processing model of 
GC formation and apply it to our sample of cluster and field dwarfs. 
In a scenario where: a) GCs are the evolved population of young 
star clusters and 6) the fraction of stars born in clusters compared to 
field stars scales with either the star formation density or gas den¬ 
sity, our results indicate that dwarfs living in galaxy clusters today 
should have an enhanced number of globular clusters per unit stel¬ 
lar mass compared to galaxies in the field with similar stellar mass, 
in agreement with observations. After accounting for stripping, we 
find that the difference in the mean specific frequency of GCs, Tjv, 
is a factor ~ 1.6, with extreme cases having Tjv as much as 10 
times larger than the field dwarf population. 

Our results lend support to the notion that dEs in clusters share 
a similar origin with isolated dirr galaxies in the field, but their 
properties have been fundamentally altered due to the cluster envi¬ 
ronment. This idea, although perhaps the most popular scenario for 
the formation of dEs, appeared to be in tension with a larger spe¬ 
cific frequency of GCs observed for cluster dwarfs. Here we find 
that such a trend could arise from the different star formation his¬ 
tories of the samples. 

We emphasize, however, that the progenitors of cluster and 
field dwarfs of similar mass at 2 = 0 can differ substantially in 
their masses at earlier times. Observational studies that are seeking 
to find a population of isolated dwarfs with the same progenitors as 
dwarfs in clusters today, may need to focus on present-day gas-rich 
galaxies that are ~ 3 times more massive than the target dE popu¬ 
lation. This difference corresponds to the freeze-out of star forma¬ 
tion that occurs for cluster dwarfs after infall, while the analog field 
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population can continue to form stars at a normal rate outgrowing 
their cluster peers from 2 ~ 1 until now. 

The results presented here have a number of predictions that 
can be tested observationally. For example, in our simulations a sig¬ 
nificant number of cluster dwarfs show intense starburst episodes, 
the majority of which are then followed by a rapid shut-off of star 
formation. Such events will leave a clear imprint in the stellar popu¬ 
lations of cluster dwarfs, likely with a-element enhancements com¬ 
pared to dwarfs in the field where star formation proceeded more 
slowly. Such a trend could be seen in observations when detailed 
metallicity data becomes available for dEs in nearby clusters. 

Another interesting prediction pertains to the kinematics of 
GCs around dEs. Our modeling of GCs suggests that some of the 
native GCs in dwarfs may have been stripped (or are in the process 
of being stripped) by the host potential. These stripped GCs will 
retain information about their origin imprinted in their phase-space 
properties (positions and velocities) that could be used to design 
observational strategies to identify the full population of native GCs 
born in dE galaxies. 

Finally, current and upcoming observations provide valuable 
insight into the structural properties of cluster dwarfs. Results so far 
have revealed a fascinating variety, with dEs showing little or no ro¬ 
tation, well defined but faded disks, and even counter-rotating com¬ 
ponents (Toloba et al 2014a,b). Moreover, dEs in the Local Group 
are also starting to show interesting differences in their star for¬ 
mation histories that could be linked to their different orbits (Geha 
et al. 2015). In a follow-up study, we will focus on the analysis 
of zoom-in simulations of clusters that have the resolution needed 
to address the morphological transformation of dwarfs in clusters. 
This will provide a comprehensive theoretical view of the evolu¬ 
tionary link between dEs and their suspected dirr progenitors and 
the many paths that have helped to shape this dichotomy in the 
dwarf galaxy population. 
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